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Abstract 



We examine signals for sparticle production at the Tevatron within the 
framework of gauge mediated supersymmetry breaking models for four dif- 
ferent model lines, each of which leads to qualitatively different signatures. 
We identify cuts to enhance the signal above Standard Model backgrounds, 
and use ISAJET to evaluate the SUSY reach of experiments at the Fermilab 
Main Injector and at the proposed TeV33. For the model lines that we have 
examined, we find that the reach is at least as large, and frequently larger, 
than in the mSUGRA framework. For two of these model lines, we find that 
the ability to identify 6-quarks and r-leptons with high efficiency and purity 
is essential for the detection of the signal. 
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I. INTRODUCTION 



Models where gauge interactions |TJ,0] rather than gravity serve as messengers of super- 
symmetry breaking have been the focus of many recent phenomenological analyses 0-[15| of 



supersymmetry (SUSY). In these models, sparticle masses and decay patterns differ from 
those in the extensively studied mSUGRA model |L6| which has served as the framework for 



many experimental analyses of supersymmetry. Perhaps the most important difference be- 
tween the mSUGRA framework and gauge-mediated SUSY breaking (GMSB) models with 
a low SUSY breaking scale is the identity of the lightest SUSY particle (LSP). In the former 
case, the lightest neutralino (Z\) is almost always the LSP, while in the GMSB framework, 
the gravitino is much lighter than other sparticles. Moreover, while the gravitino is essen- 
tially decoupled in mSUGRA scenarios, the couplings of the Goldstino (which forms the 
longitudinal components of the massive gravitino), though much smaller than gauge cou- 
plings, may nonetheless be relevant for collider physics in that they can cause the next to 
lightest SUSY particle (NLSP) to decay into a gravitino inside the detector. The precise 
decay pattern and lifetime of the NLSP depends on its identity and on model parameters. 
For instance, if Z\ is the NLSP, it would decay via Z\ — ■> G + 7, and if kinematically allowed, 
also via Z\ — > G + Z, or into the various Higgs bosons of SUSY models via Z\ — > G + h, H, 
or A. If, on the other hand, the NLSP is a slepton, it would decay via i — > G + £, etc. 
Sparticles other than the NLSP decay only very rarely to gravitinos, so that it is safe to 
neglect these decays in any analysis. Thus heavier sparticles cascade decay as usual to the 
NLSP, which then decays into the gravitino as described above. 

Sparticle signatures differ from those in the mSUGRA framework for two basic reasons. 
First, if the NLSP is not the lightest neutralino, the cascade decay patterns to the NLSP 
are modified. Second, the NLSP (which need not be electrically neutral) itself decays into 
a gravitino and Standard Model (SM) particles. The gravitinos escape the experimental 
apparatus undetected resulting in J£t in SUSY events. In the GMSB framework, however, 
neutralino NLSP decays may also result in isolated photons or Z or Higgs bosons which could 
provide additional handles to reduce SM backgrounds to the SUSY signal. If the NLSP is a 
slepton, all SUSY events should contain leptons of the same flavour as the slepton NLSP in 
addition to Ifo. While it may be possible to have other candidates for the NLSP, this does 
not seem to be the case in the simplest realizations of the GMSB framework, and we will 
not consider this possibility any further. 

Within the minimal GMSB framework, supersymmetry breaking in a hidden sector is 
communicated to the observable sector via SM gauge interactions of messenger particles 
(with quantum numbers of SU{2) doublet quarks and leptons) whose mass scale is charac- 
terized by M. As a result, the soft SUSY breaking masses induced for the various sparticles 
are directly proportional to the strength of their gauge interactions. Thus, coloured squarks 
are heavier than sleptons, and gluinos are heavier than electroweak gauginos. The observ- 
able sector sparticle masses and couplings are determined (at the scale M) by the GMSB 
model parameter set, 

A, M, n 5 , tan/?, sign(fi), C grav . (1.1) 

Of these, A is the most important parameter in that it sets the scale of sparticle masses. 
The model predictions for the mass parameters at the scale M are then evolved down to 



2 



the sparticle mass scale via renormalization group evolution (RGE). Radiative breaking of 
electroweak symmetry determines \fi\. The weak scale SUSY parameters depend only weakly 
on the messenger mass scale M, since this primarily enters as the scale at which the mass 
relations predicted by the model are assumed to be valid. There is an additional dependence 
of the sparticle spectrum on M due to threshold effects jr7||, but this is also weak as long 
as M/A is not very close to unity. Messenger quarks and lepton, it is assumed, can be 
classified into complete vector representations of SU(5): the number (715) of such multiplets 
is required to be < 4 for the messenger scale M = (9(100 TeV) in order that the gauge 
couplings remain perturbative up to the grand unification scale. Finally, the parameter 
Cgrav > 1 Jl3| (essentially, the ratio of hidden sector to messenger sector SUSY breaking 
vevs) can be used to dial the gravitino mass beyond its minimum value. Effectively, C grav 
parametrizes the rate for sparticle decays into a gravitino. This decay is most rapid when 
Cgrav = 1, while for larger values of C grav the NLSP may decay with an observable decay 
vertex, or may even be sufficiently long-lived to pass all the way through the detector. In 
this extreme case, SUSY event topologies would be identical to those in the mSUGRA model 
if Z\ is the NLSR However, for the case where a charged slepton is the NLSP, SUSY events 
would necessarily contain a pair of penetrating tracks from the long-lived slepton NLSP, 
which might be detectable at the Tevatron as "additional (possibly slow) muons" ]TT|,P 



In our analysis, we assume that the NLSP decays promptly and fix C grav = 1; i.e. we do 
not attempt to model the additional handle displaced vertices might provide to reduce SM 
backgrounds. 

Many of the phenomenological implications depend only weakly on the parameter M. 
Thus, the A — tan (3 plane provides a convenient panorama for illustrating the diversity of 
phenomenological possibilities in GMSB scenarios. This is shown in Fig. 1 for M = 3A 
and a) n 5 = 1, b) n 5 = 2, c) n 5 = 3 and d) n 5 = 4. We choose ji to be positive since 
for this choice the model predictions |6|JT8[| are well within experimental constraints [l9j 
from the decay b —>■ s'j over essentially the whole plane. In Region 1 in Fig. la-d (the 
boundaries of these regions are the heavy solid lines), the lightest neutralino is the NLSP, 
so that Z\ — > G'j (and to Z and Higgs bosons if these decays are kinematically allowed). 
In Region 2, < , while other sleptons are heavier than Zi, and cascade decays of 
sparticles terminate in f±, except immediately above the boundary between Regions 1 and 
2 where the decay Z\ — > f±r is kinematically forbidden. For parameters in Region 2, we 
thus expect an excess of r leptons in SUSY events ||. In Regions 3 and 4 in Fig. lb-d, not 
only f\, but also e.\ ~ e~R and jli ~ fl^, are lighter than Z\. Thus neutralinos are effectively 
sources of real (dominantly right-handed) sleptons. In Region 3, l\ — > £G because its decay 
£i — > f\rl (£ = e, fi) is kinematically forbidden. The decay /ti — > u^fiu T which occurs 
via suppressed muon Yukawa couplings is kinematically allowed, and may compete with 
the decay to jx\ — > p,G; for C grav = 1, we find that this decay (which has been included 
in our computation), is unimportant. In Region 4, the decays l\ — > f\fl and £\ — > f\rl 
are also allowed, and compete with the gravitino decay of i\. Frequently, the stau decays 
of i\ dominate its decays to gravitino, and then, as for Region 2, SUSY events will be 
characterized by an abundance of taus in the final state. Signals from sparticle production 
will clearly depend on which of these regions the model parameters happen to lie. 

The grey regions in Fig. 1 are excluded because the observed pattern of electroweak 
symmetry breaking is not obtained: in the wedge in the upper left corner, m\ R < while 
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in the band on top, m\ < 0. The non-observation of sparticle signatures in experiments at 
LEP excludes other portions of the plane. Within the MSSM charginos have been excluded 
if — 90 — 95 GeV. While this limit has been obtained assuming that charginos and 
selectrons decay into a stable neutralino which escapes detection, we expect that an even 
more striking signature is obtained if Z\ decays via Z\ — ► jG. The leftmost dot-dashed line 
in Fig. 1 is the contour = 95 GeV: to its left, charginos are lighter than 95 GeV. To 
assist the reader in assessing the sparticle mass scale, we have also shown mass contours 
for = 200 GeV and = 350 GeV. MSSM searches for acollinear electron pairs 

exclude selectrons lighter than about 90 GeV. This limit should certainly be valid within 
this framework if e,\ — ► Ge (Region 3) or even if it decays to Z 1 that subsequently decays to 
a photon (Region 1). The dotted line is the contour = 90 GeV. For the case where e\ 
mainly decays to fi, the actual bound may be somewhat weaker, and closer to the MSSM 
stau bound ~ 76 GeV. The most stringent experimental limit for the n 5 = 1 and n 5 = 2 
cases in frames a) and b) comes from the the LEP search [^U| for 77+ fir events from 
e + e~ — > ZxZi production. The cross section for this process depends on the selectron mass. 
The ALEPH analysis for 715 = 1 results in the lower limit, m-^ > 84 GeV ||22|| . For 

larger values of n 5 , the selectron to Z\ mass ratio is smaller, so that the corresponding cross 
section is even larger than in the 77,5 = 1 case. Indeed the DELPHI collaboration |20| has 
obtained a preliminary bound <i 88 GeV for = 2, for parameters in Region 1. If 

m fl < , Zis act as sources of staus and add to the signal from direct stau pair production. 
The DELPHI search for acollinear tau pairs still limits < 86 — 90 GeV, and also the 
bounds m fl < 76 GeV, regardless of m^. In Fig. 1, in the horizontally hatched region 
m z 1 — 84 GeV, whereas in the region with vertical hatches, < 76 GeV. Finally, the 
LEP experiments [24] have a preliminary bound of 91-95 GeV on the mass of the SM Higgs 
boson. Since, for small values of tan j3 the lighter Higgs boson h of the GMSB framework 
is frequently close to the SM Higgs boson, we also show the regions with rrih < 95 GeV 
(the diagonally hatched area in the lower left corner) in Fig. 1. Furthermore, LEP analyses 
exclude ttia ~ 83 GeV when tan (3 is large. This excludes the thin (diagonally hatched) sliver 
where tan (3 ~ 53. The reader should appreciate that the various shaded regions that we 
have shown are not formal experimental limits, but indicate the reach of present experiments 
within the GMSB framework. 

We see from Fig. 1 that current experiments have already probed Regions 1 and 2 if 
n.5 > 2. On the other hand, for n.5 = 1, we have just these two regions, while n§ = 2, all four 
regions are still possible. Since experimental signatures within the GMSB framework differ 
significantly from those in the MSSM and mSUGRA models, it is of interest to reassess the 
sensitivity of Tevatron experiments to signals from sparticle production at the upcoming 
Run II of the Tevatron Main Injector (MI) as well as at the proposed luminosity upgrade 
(dubbed TeV33) where an integrated luminosity ~ 25 fb~ x might be accumulated. This is 
the main purpose of this paper. We had begun this program in an earlier study where 
we had computed cross sections for various SUSY event topologies for models with n 5 = 1 
expected at the Tevatron: in this case, the NLSP is dominantly the hypercharge gaugino. 
Here, we first repeat this analysis for somewhat different model parameters, using cuts and 
acceptances more appropriate to Run II. We have also fixed a bug in the program |26[ which 
resulted in an underestimate of the chargino pair production cross section. We also examine 
cases with larger values of for which we expect the phenomenology to change qualitatively 
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from our earlier study. Toward this end, we first examine a model line with = 2 with 
tan/3 = 15 where fx is the NLSP and significantly lighter than other sleptons. Next, we 
examine a model line with n 5 = 3 where all right handed sleptons are roughly degenerate in 
mass (the co-NLSP scenario), and where e x (/ti) dominantly decay via e.\ — > eG (fit — * (jlG). 
Finally, we examine a non-minimal model where the NLSP is dominantly a Higgsino-like 
neutralino. This is not because we believe this is any more likely than the mGMSB scenarios 
previously discussed, but because it leads to qualitatively different experimental signatures. 
In view of the fact that the underlying mechanism of SUSY breaking, and hence the re- 
sulting mass pattern, is unknown it seems worthwhile to explore implications of unorthodox 
scenarios, particularly when they lead to qualitative differences in the phenomenology. 

In the next Section, we describe the upgrades that we have made to ISAJET to facilitate 
the simulation of the minimal GMSB framework that we have described, as well as several 
of its non-mimimal extensions. In Sec. Ill we specify four different model lines and discuss 
strategies for separating the SUSY signal from SM background for each of these. Our main 
result is the projection for the reach of experiments at the MI and at TeV33. We end in 
Sec. IV with a summary of our results together with some general remarks. 



II. SIMULATION OF GAUGE MEDIATED SUSY BREAKING SCENARIOS 

We use the event generator program ISAJET v 7.40 for simulating SUSY events at the 



Tevatron. Since ISAJET has been described elsewhere |]2"5" |, we will only discuss recent im- 
provements that we have made that facilitate the simulation of the mGMSB model specified 
by the parameter set (|1 . 1| ) , and also, some of its variants. The 'GMSB option' allows one 
to use the parameter set (|1.1|) as an input. ISAJET then computes sparticle masses at the 
messenger scale M, then evolves these down to the lower scale relevant for phenomenology, 
and finally calculates the 'MSSM parameters' that are then used in the evaluation of spar- 
ticle cross sections and decay widths. The decays of neutralinos into gravitinos, Zi — > Gj, 
Zi — > GZ and Zi — > Gh, H, A as well as (approximately) the Dalitz decay Zi — > e + e~G 
are included in ISAJET. The decays t\ — > £G and f\ — > tG, as well as the three body 



decays [^j l\ — > t\tI and i\ — > fir£, which are mediated by a virtual neutralino, have 
also been included. The widths for corresponding three body decays mediated by virtual 
chargino exchange are suppressed by the lepton Yukawa coupling, and are also included. 
These decays can be significant only for smuons, and only when — m fl < m T so that the 
neutralino-mediated three body decays of fx\ are kinematically very suppressed or forbidden. 
Although for C grav = 1 we have not found this to be important, for larger values of C grav , 
the (long-lived) smuon may dominantly decay via the chargino mediated decay to a stau, 
and may alter the apparent curvature of the 'smuon track' in the detector ||28|| . 

We have also included in ISAJET the facility to simulate several non-minimal gauge- 
mediated SUSY breaking models that involve additional parameters. While these will be 
irrelevant to the analysis performed in our present study, we have chosen to describe this 
for completeness because it may prove useful to readers studying extensions of the minimal 
class of models. 



The parameter fl allows the user to adjust [|3D| the ratio between the gaugino and 
scalar masses by scaling the former by the factor fl which is equal to unity in the 
minimal GMSB framework. 
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In GMSB models, additional interactions are needed to generate the dimensional fi 
and B parameters that are essential from phenomenological considerations. These 
interactions can split the soft SUSY breaking masses of Higgs and lepton doublets 
(at the messenger scale) even though these have the same gauge quantum numbers. 
These additional contributions to the squared masses of Higgs doublets that couple to 
up and down type fermions, are parametrized by 5m 2 Hu and 5m 2 H , respectively. 
These parameters are zero in the minimal model. 

If the hypercharge D-term has a non-zero expectation value Dy in the messenger 
sector, it will lead to additional contributions to sfermion masses at the messenger 
scale which may be parametrized |3IJ as 5m 2 j = g'YDy. The value of Dy (which is 
zero in the minimal GMSB framework) is constrained as it can lead to an unacceptable 
pattern of electroweak symmetry. 



• Finally, allowing incomplete messenger representations |T7| can effectively result in 
different numbers (n 5l , n 52 and n 5s ) for each factor of the gauge group. 

ISAJET allows the user to simulate these non-minimal models using the GMSB2 command. 

To model the experimental conditions at the Tevatron, we use the toy calorimeter simu- 
lation package ISAPLT. We simulate calorimetry covering —4 < 77 < 4 with a cell size given 
by At] x A0 = 0.1 x 0.087, and take the hadronic (electromagnetic) calorimeter resolution 
to be 0.7 /VE (0.15 /\/~E). Jets are defined as hadronic clusters with E T > 15 GeV within a 
cone of AR = Ar/ 2 + Acp 2 = 0.7 with \r]j\ < 3.5. Muons and electrons with Et > 7 GeV 
and 1 77^ I < 2.5 are considered to be isolated if the the scalar sum of electromagnetic and 
hadronic Ej- (not including the lepton, of course) in a cone with AR = 0.4 about the lepton 
to be smaller than max(2 GeV, Er(£)/4:). Isolated leptons are also required to be separated 
from one another by AR > 0.3. We identify photons within |r/ 7 | < 1 if Et > 15 GeV, 
and consider them to be isolated if the additional Et within a cone of AR = 0.3 about the 
photon is less than 4 GeV. Tau leptons are identified as narrow jets with just one or three 
charged prongs with pt > 2 GeV within 10° of the jet axis and no other charged tracks in 
a 30° cone about this axis. The invariant mass of these tracks is required to be < m T and 
the net charge of the three prongs required to be ±1. QCD jets with Et = 15(> 50) GeV 
are misidentified as taus with a probability of 0.5% (0.1%) with a linear interpolation in 
between. Finally, for SVX tagged 6-jets, we require a jet (satisfying the above jet criteria) 
within \r)j\ < 1 to contain a .B-hadron with pr > 15 GeV. The jet is tagged as a 6-jet with 
a probability of 55%. Charm jets (light quark or gluon jets) are mistagged as 6-jets with a 
probability of 5% (0.2%). 



III. THE REACH OF TEVATRON UPGRADES FOR VARIOUS MODEL LINES 

Within the GMSB framework, sparticle signatures, and hence the reach of experimental 
facilities, are qualitatively dependent on the nature of the NLSP. Here, we examine the reach 
of experiments at the Tevatron Main Injector as well as that of the proposed TeV33 upgrade 
for four different model lines [31] where the NLSP is (A) dominantly a hypercharge gaugino, 



(B) the stau lepton, fj., with other sleptons significantly heavier than f±, (C) again the stau, 
but ex and fix are essentially degenerate with fx, and (D) dominantly a Higgsino. We fix the 



6 



messenger scale M = 3A, fi > and C grav = 1 throughout our analysis. We use ISAJET to 
compute signal cross sections, incorporating cuts and triggers to simulate the experimental 
conditions at the Tevatron together with additional cuts that serve to separate the SUSY 
signal from SM backgrounds. We project the reach of future Tevatron upgrades for each of 
these scenarios. 



A. Model Line A: The Bino NLSP Scenario 

We see from Fig. la that for n 5 = 1, the lightest neutralino is the NLSP as long as tan/3 
is not very large. Since the value of computed from radiative breaking of electroweak 
symmetry is rather large, the NLSP is mainly a bino. To realize the bino NLSP model line, 
we fix tan/5 = 2.5 which ensures that sleptons are significantly heavier than m^^. Sparticles 

cascade decay to Z\ which then mainly decays via Z\ — > 7G. Thus almost all SUSY events 
contain at least two hard isolated photons. 

In Fig. we show the mass spectrum of sparticles that might be in the Tevatron range 
versus A, which sets the sparticle mass scale, while in frame b) we show the cross sections for 
the most important sparticle production mechanisms at the Tevatron. We see that chargino 
pair production and W\Z 2 production dominate because squarks and gluinos are beyond 
the Tevatron reach. The production of right-handed slepton pairs is suppressed relative 
to chargino/neutralino production by over an order of magnitude. Values of A smaller 
than ~ 70 TeV are excluded by the LEP search for 77+ fir events. For A < 80 TeV 
(corresponding to < 90 — 100 GeV), the two body decay Wi — > WZi is kinematically 

suppressed, and the chargino mainly decays via W\ — > T\v T or W\ — ► qqZi] for A > 80 TeV, 
the decay W\ — > W Z\ dominates. The neutralino Z 2 dominantly decays via Z 2 — > Zih 
(for A > 90 TeV) when this decay is not kinematically suppressed: otherwise it decays 
via Z 2 — > £i£, with roughly equal branching fractions for all three lepton flavours. We 
thus expect that W1W1 and Z 2 W\ production will mainly lead to jetty events (counting 
hadronically decaying taus as jets) possibly with additional e and \i plus photons plus Et- 
We use ISAJET to classify the supersymmetric signal events primarily by the number of 
isolated photons — events with < 2 photons arise when one or more of the photons is outside 
the geometric acceptance, has too low an Et, or happens to be close to hadrons and thus 
fails the isolation requirement. We further separate them into clean and jetty events and 
then classify them by the number of isolated leptons (e and /1). In addition to the acceptance 
cuts described in Sec. 2, we impose an additional global requirement > 40 GeV, which 
together with the presence of jets, leptons or photons may also serve as a trigger for these 
events. 

Before proceeding to present results of our computation, we pause to consider SM back- 
grounds to these events. We expect that the backgrounds are smallest in the two photon 
channel, which we will mainly focus on for the purpose of assessing the reach. We have not 
attempted to assess the background because the recent analysis by the DO collaboration PTJ . 
searching for charginos and neutralinos in the GMSB framework, points out that the major 
portion of the background arises from mismeasurement of QCD jets and for yet higher values 
of Et from misidentification of jets/leptons as photons. In other words, this background 
is largely instrumental, and hence rather detector-dependent. From Fig. 1 of Ref. [21], we 
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estimate the inclusive 27+ Et > 40 GeV (60 GeV) background level (for £^(71,72) > 
(20 GeV, 12 GeV)) to correspond to ~ 0.9 (0.1) event in their data sample of ~ 100 pb^ 1 . 
The background from jet mismeasurement, of course, falls steeply with Et- The inclusive 
27+ Et background is also sensitive to the minimum Et of the photon. 

To assess how changing the photon and Et requirements alter the SUSY signal, in Fig. |3| 
we show the signal distribution of (a) £^(72), the transverse energy of the softer photon in 
two photon events, and (b) Et in 77+ Et events that pass our cuts, for three values of A. 
The following is worth noting. 

• For A ~ 100 TeV (which we will see is in the range of the Tevatron bound), reducing 
the Et(j) cut does not increase the signal. In fact, it may be possible to further 
harden this cut to reduce the residual backgrounds. Although we have not shown it 
here, we have checked that increasing the cut on the hard photon to £t(7i) > 40 GeV 
results in very little loss of signal for A > 100 TeV. 

• In view of our discussion about SM backgrounds, it is clear that requiring Et > 60 GeV 
greatly reduces the background with modest loss of signal. Indeed, it may be possible 
to reduce the background to negligible levels by optimizing the cuts on E T of the 
photons and on Et- 

The results of our computation of various topological cross sections at a 2 TeV pp collider 
after cuts are shown in Fig. |] for (a) photon, (b) one photon, and (c) two photon events. 
In this figure, we have required that Et > 60 GeV. As mentioned, this reduces the cross 
section by just a small amount, especially for the larger values of A in this figure. The solid 
lines correspond to cross sections for events with at least one jet, while the dashed lines 
correspond to those for events free of jet activity. The numbers on the lines denote the 
lepton multiplicity, and are placed at those A values that we explicitly scanned. Finally, 
the heavy solid line represents the sum of all the topologies, i.e. the inclusive SUSY cross 
section after the cuts. We note the following. 

• We have comparable signal cross sections in I7 and 27 channels. Since the background 
in the latter is considerably smaller (recall a significant portion of it is from fake 
photons), the maximum reach is obtained in the 27 channel. 

• As anticipated, events with at least one jet dominate clean events, irrespective of the 
number of photons. 

We may obtain a conservative estimate of the reach by assuming an inclusive 27+ Et > 
60 GeV background level of 0.1 event per 100 pb' 1 ; i.e. assuming a background level of 

1 fb. This corresponds to a "5cr reach" of 3.5 fb (1 fb) for an integrated luminosity of 

2 fb- 1 (25 fb- 1 ) at the Tevatron, or A < 110 TeV (130 TeV) at the Main Injector (TeV33 
upgrade). As we have mentioned, it may be possible to further reduce the background by 
hardening the £r(7) arid Et requirements with only modest loss of signal. The background 
may also be reduced if jet/lepton misidentification as a photon is considerably smaller than 
in Run I ETJ. If we optimistically assume that the reach is given by the 5 (10) event level 
at the Main Injector (TeV33), we would be led to conclude that experiments may probe 
A values as high as 118 TeV (145 TeV) at these facilities. It should be remembered that 
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A = 118 TeV corresponds to m~ g ~ 950 GeV, almost equal to what is generally accepted as 
the qualitative upper limit from fine tuning arguments. 

We see from Fig. that the l\£\ production cross section exceeds 1 fb for A < 100 TeV. 
Since slepton production can lead to spectacular ££77+ E T events of the type observed by 
the CDF Collaboration ||23|| , it appears reasonable to ask whether signal from slepton pair 
production might be observable at TeV33, and further, whether it can be separated from 
a similar signal from chargino pair production when each W\ — > £vZ\ — ► IvG^ |32|| . The 
SM physics backgrounds come from WW'yy production which for E^ > 10 GeV has a 
production cross section of 0.15 ± 0.05 fb, so that ~ 0.1 such event is expected in a 
data sample of 25 fb -1 . The background from tt production is estimated to be even smaller. 
In Fig. [| we show the total cross section for clean ££77+ Et events after cuts (solid) and 
the corresponding cross section from just e.\ and \i\ pair production (dashed). We see that 
for A < 115 TeV (corresponding to ~ 200 GeV), five or more signal events should 
be present at TeV33, with about 60% of these having their origin in direct production of 
sleptons. Slepton pair production alone yields five events for < 180 GeV. If instrumental 
backgrounds from jets faking an electron or photon turn out to be negligible, direct detection 
of sleptons as heavy as 180 GeV may be possible at TeV33 p3[ for model line A. 



B. Model Line B: The Stau NLSP Scenario 

From Fig. 1, we see that we can obtain fi as the NLSP for a wide range of GMSB 
parameters. Here, we choose = 2, and take tan/? = 15 to make e\ and fi\ somewhat 
heavier than f\, with other parameters as before. In Fig. ^] we show a) relevant sparticle 
masses, and b) cross sections for the main sparticle production mechanisms versus A. For 
A > 30 TeV, mfj < but for A < 32 TeV, Z\ — > rf\ is kinematically forbidden, and Z\ 

would decay via the four body decay Z\ — > v T fiW* (which is not yet included in ISA JET) 
or via its photon mode considered above. In our study, we only consider A > 35 TeV, the 
region safe from LEP constraints. Gluinos and squarks are then too heavy to be produced 
at the Tevatron, and sparticle production is dominated by chargino, neutralino and, to a 
lesser extent, slepton pair production. 

The two body decay W\ — > f\V is always accessible, while the decay W\ — > WZ\ becomes 
significant only for A > 45 TeV ^ 210 GeV). The branching fraction for Z 2 decays 

are shown in Fig. ^j\a. We see that Z% decays via f\T with a branching fraction that exceeds 
0.5 if m^ 2 < 300 GeV, and further that branching fractions for Z 2 — > l\l (£ = e, /x) are not 

negligible. For the value of tan/3 in this Figure, the decay Z 2 — > Z\h is only important for 
relatively large values of A. The lightest neutralino Z\ mainly decays via Z\ — > t\t, though 
for large enough values of A its decay to sleptons of other families may also be significant. 
The decay pattern of the lighter selectron and smuon are illustrated in Fig. [7|&. For small 
values of A (Region 2 of Fig. 1), the decay l\ —>■ IZ\ is kinematically allowed and dominates. 
For larger values of A (Region 4 of Fig. 1), where this channel is closed, the neutralino is 
virtual and £1 — > t^^t^ or t[ — > t{1~t + . These decays dominate the decay i\ — > £G. 
The upshot of these decay patterns is that SUSY events may contain several tau leptons 
from sparticle cascade decays. At the very least, each event will contain a pair of rs (in 
addition to other leptons, jets and E T ) from fx produced at the end of the decay cascade. 
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It is worthwhile to note that the two rs could easily have the same sign of electric charge. 

The observability of SUSY realized as in this scenario thus depends on the capability of 
experiments to identify hadronically decaying tau leptons, and further, to distinguish these 
from QCD jets. Following the same logic as in the = 1 case above, we now classify SUSY 
events by the number of identified taus, and further separate them into jetty and clean event 
topologies labelled by the number of isolated leptons (e and //). It should be remembered 
that the efficiency for identifying taus is expected to be smaller than for identifying photons 
- first, the tau has to decay hadronically, and then the hadronic decay products have to 
form a jet. 

In our analysis, in addition to the basic acceptance cuts discussed in Section II, we require 
f*T > 30 GeV together with at least one of the following which serve as a trigger for the 
events: 

• one lepton with pr(£) > 20 GeV, 

• two leptons each with pt{£) > 10 GeV, 

• $ T > 35 GeV. 

In addition, we also impose the additional requirements: 

• a veto on opposite sign, same flavor dilepton events with Mz — 10 GeV < m(££) < 
Mz + 10 GeV to remove backgrounds from WZ and ZZ and high Z production, 
and 

• for dilepton events, require A<ft(££') < 150° (£, £' = e, r) to remove backgrounds from 
Z — > rf events. 

The dominant physics sources of SM backgrounds to n-jet + m-leptons + $t events, 
possibly containing additional taus, are W, 7* or Z + jet production, tt production and 
vector boson pair production. Instrumental backgrounds that we have attempted to estimate 
are $t from mismeasurement of jet energy and mis-identification of QCD jets as taus. 

We have checked that even after these cuts and triggers, SM backgrounds from W pro- 
duction swamp the signal in channels with no leptons or just one identified lepton (e, /i or 
r). The former is the canonical J£t signal, which after optimizing cuts, may be observable 
at Run 2 if gluinos are lighter than ~ 400 GeV. We do not expect that this signal from 
gluino and squark production will be detectable since even for A = 35 TeV, irtg = 578 GeV 
with squarks somewhat heavier. For this reason, and because there are large single lepton 
backgrounds from W production, we focus on signals with two or more leptons in our study. 
Also, because the presence of r's is the hallmark of this scenario, we mostly concentrate on 
leptonic events with at least one identified r. 

We begin by considering the signal and background cross sections for clean events. These 
are shown in Table |. Events are classified first by the number of identified taus, and then 
by the lepton multiplicity; the C in the topology column denotes "clean" events. For each 
topology, the first row of numbers denotes the cross sections after the basic acceptance 
cuts and trigger requirements along with the Z veto and the A<f) cut discussed above. We 
see that there is still a substantial background in several of the multilepton channels. This 
background can be strongly suppressed, with modest loss of signal by imposing an additional 
requirement, 
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• PTvis{ri) > 40 GeV, 

on the visible energy of the hardest tau in events with at least one identified tau. In the 
background, the rs typically come from vector boson decays, while in the signal a substantial 
fraction of these come from the direct decays of charginos and neutralinos that are substan- 
tially heavier than Mz (remember even = 103 (132) GeV for A = 40 (50) TeV): thus 
signal taus pass this cut more easily. A few points about this Table are worth mentioning. 

1. The signal cross sections in each channel are at most a few fb, and with an integrated 
luminosity of 2 fb' 1 , the individual signals are below the 5a level even for A = 40 TeV. 
For the integrated luminosity expected in Run II of the MI we will be forced to add 
the signal in various channels and see if this inclusive signal is observable. 

2. The sum of the signal in all the channels in Table [|, except the 1t1£ channel which 
has a very large background, is shown in the next two rows both with and without 
the pt cut on the r, while in the last two rows we list a(sig) / \Ja(back). We see that 
a somewhat better significance is obtained after the Ptws(ti) > 40 GeV cut. 

3. We see that the inclusive SUSY signal in the clean channels for the A = 40 TeV case 
should be detectable with the Run II integrated luminosity, whereas for the A = 50 TeV 
case an integrated luminosity of 12 fb -1 is needed for a 5cr signal. 

4. We caution the reader that about 25-30% of the r background comes from mis-tagging 
QCD jets as taus (except, of course, for the W backgrounds and the backgrounds in the 
C2t£ channels which are almost exclusively from these fake taus). Thus our estimate 
of the background level is somewhat sensitive to the r faking algorithm we have used. 
The signal, on the other hand, almost always contains only real rs, so that improving 
the discrimination between r and QCD jets will lead to an increase in the projected 
reach of these experiments. 

5. In some channels the background is completely dominated by fake taus. For instance, 
after the prvisij) cut, the C\r\l background from W sources of just real taus is only 
1.9 fb, while the signal and other backgrounds remain essentially unaltered from the 
cross sections in Table |. Thus if fake r backgrounds can be greatly reduced, it may 
be possible to see a signal via other channels. 

Next, we turn to jetty signals for Model Line B. Cross sections for selected signal topolo- 
gies together with SM backgrounds after the Pt{ t i) — 40 GeV cut are shown in Table [TI|. 
The other topologies appear to suffer from large SM backgrounds and we have not included 
them here. 

The following features are worth nothing. 

1. We see that Model Line B results in smaller cross sections in jetty channels. This 
should not be surprising since electroweak production of charginos, neutralinos and 
sleptons are the dominant SUSY processes, and because staus are light, branching 
fractions for hadronic decays of W\ and Zi j2 tend to be suppressed. 
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2. We see from Table |I| that with the present set of cuts, not only is the signal below 
the level of observability in any one of the channels, but also that the inclusive signal 
is not expected to be observable at the MI even for the A = 40 TeV case. With an 
integrated luminosity of 25 fb^ 1 the signal for the A = 50 TeV case is observable at 
the 6<T level. 

3. As for the clean lepton case, a significant portion of the background comes from QCD 
jets faking a tau. The fraction of events with a fake tau varies from channel to channel, 
but for the 2t1£ channel in Table almost 60% of the background (in contrast to 
essentially none of the signal) involves at least one fake r. 

4. A major background to the jetty signal comes from it production. To see if we could 
enhance the signal relative to this background we tried to impose additional cuts to 
selectively reduce the top background. Since top events are expected to contain hard 
jets, we first tried to require E T (j) < 50 GeV. We also, independently, tried vetoing 
events where the invariant mass of all jets exceeded 70 GeV. While both attempts lead 
to an improvement of the signal to background ratio, the statistical significance of the 
signal is not improved (and is even degraded) by these additional cuts. We do not 
present details about this for the sake of brevity. It may be possible to reduce the top 
background by vetoing events with tagged 6-jets, but we have not attempted to do so 
here. 

For Model Line B, it appears that experiments at the MI should be able to probe A 
values up to just beyond 40 TeV in the inclusive clean multilepton channels. It appears, 
however, that it will be essential to sum up several channels to obtain a signal at the 5<r 
level. Confirmatory signals in inclusive jetty channels may be observable at the 3.7a level. 
Of course, for an integrated luminosity of 25 fb^ 1 it may be possible to probe A = 50 TeV 
even in the unfavoured jetty channels, and somewhat beyond in the clean channels. The 
situation is summarized in Fig. || where we show the signal cross sections summed over 
the selected channels for events without jets (dashed) and for events with jets (solid). The 
horizontal lines denote the minimum cross section needed for the signal to be observable 
at the 5cx level, for the two assumptions about the integrated luminosity. We note that, 
in some channels, a susbtantial fraction of background events come from QCD jets faking 
a tau — our assessment of the Run II reach is thus sensitive to our modelling of this jet 
mis-tag rate. By the same token, if this rate can be reduced, the reach may be somewhat 
increased. Finally, we remark that even though it appears that the range of parameters 
that is accessible to experiments at the MI is very limited (A < 42 TeV), these parameters 
correspond to charginos as heavy as 192 GeV and gluinos and squarks around 700 GeV. 

C. Model Line C: The co-NLSP Scenario 

This scenario can simply be obtained by choosing parameters in Regions 3 of Fig. 1, 
so that ei, fix and f\ are all approximately degenerate, and t\ and jli cannot decay to f\. 
Here, we choose n 5 = 3 and tan j3 = 3 with other parameters as before. In Fig. |S| we show a) 
relevant sparticle masses, and b) cross sections for the main sparticle production mechanisms 
at the Tevatron versus A. Aside from the fact that lighter sleptons of all three flavours have 
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essentially the same mass, the main difference from the earlier model lines that we have 
studied is that, while charginos and neutralinos dominate for lower values of A, slepton 
pair production is the dominant production mechanism for A > 40 TeV (corresponding to 
mi> 125 GeV). 

In Fig. [H] we show the branching fractions for a) chargino, and b) neutralino decays 
versus A. For small A in frame a), charginos dominantly decay via W\ — ► f\v r since the cor- 
responding decays to smuons and selectrons are suppressed by the lepton Yukawa coupling. 
As A increases, decays to sneutrinos and the heavier (dominantly left-handed) sleptons be- 
come accessible. Since these occur via (essentially unsuppressed) gauge interactions, these 
rapidly dominate the decay to f\. The decay W\ — > Z{W also becomes significant for 
m Wi — GeV. Turning to Z 2 decays shown in frame b), we see that these dominantly 
decay to sleptons with branching fractions more or less independent of the lepton flavour. 
Again, since Z 2 is dominantly an 577(2) gaugino, decays to the heavier (dominantly left- 
handed) sleptons and sneutrinos dominate when these are kinematically unsuppressed. The 
branching fraction for the decay Z 2 — > Z\h is also significant. From the plot of sparti- 
cle masses in Fig. |9]a, we see that the heavier charged sleptons and sneutrinos decay via 
ji — ► jZ\ while fi — > fG (/ = £,v). The lightest neutralino decays via Z\ — > l\l with 
branching fractions essentially independent of the lepton flavour. 

The bottom line of these decay patterns is that even though f\ is strictly speaking the 
NLSP, we expect a large multiplicity of isolated leptons (e and \i) from sparticle production 
at colliders. This is because all flavours of sleptons are roughly equally produced in SUSY 
cascade decays, and decays of e.\ and \x\ do not involve a stau at an intermediate stage. This 
is illustrated in Fig. [0] for A = 30 TeV and A = 40 TeV, where we show the multiplicity 
distributions for both n e + n M and n e + + n T , for events satisfying the basic acceptance 
cuts (see Section II) and trigger conditions, but not the additional requirements described 
in the Section IIIB. We see that while the lepton multiplicity is largest for two leptons (due 
to production of Ir pairs), a very sizeable fraction of signal events have both n e + and 



n e + rifi + n T > 4 for which backgrounds from SM sources, shown in Table |HI| , are very small. 
Here, for the n e + + n r > 4 background sample, we found that all the background events 
that passed the cuts automatically satisfied n e + n^> 2. We also impose this requirement, 



which should facilitate triggering on these events even without a r trigger, on the signal |34 

In our simulation we found that the background from W and Z plus jet production, as 
well as tt production are negligible. We see that the 4-lepton background is an order of 
magnitude smaller than the corresponding 3£ backgrounds in Tables | and |IJ (even though 
the cuts there are more stringent than in Table p|) . In contrast, for the signal we see from 
Fig. |TI| that the rate for 3£ events is much smaller than that for > 4£ events. Thus the > 4£ 
channel offers the best hope for detection of the SUSY signal for model line C. 

The SUSY reach for the co-NLSP model line is illustrated in Fig. [12] where we show the 
signal cross section versus A for inclusive events with n e +n^ > 4 (dashed) and n e +n^+n T > 4 
(solid) where we require, in addition, n e +n^ > 2. Here we have summed the cross section for 
events with and without jets as this offers the greatest reach. The corresponding horizontal 
lines denote the levels where the signal will be just detectable at the '5cr level' (with a 
minumum of at least five events) at the MI and at the proposed TeV33 upgrade. We 
observe the following. 

• With an integrated luminosity of 2 fb -1 , the signal is rate-limited in the n e + n„ > 4 
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channel, and experiments at the MI should be able to probe A out to about 45 TeV 
(corresponding to charginos heavier than 300 GeV) if we require a minimum signal level 
of five events. Including taus increases the signal but also increases the background so 
that the reach is only marginally improved. Since this background largely comes from 
tau mis-identification, it should be kept in mind that our projection for the reach via 
the n e + n M + n T > 4 channel is somewhat dependent on our simulation of this. 

• For an integrated luminosity of 25 fb~ x we see that the projected increase in the 
background in the channel that includes taus actually leads to a reduction of the 
reach, and the greatest reach is obtained via events with n e + > 4 leptons for 
which the background is very small. In our assessment of this reach we have assumed 
that backgrounds from hadrons or jets mis-identified es and /is are negligible. The 
reach of TeV33 experiments should then extend to A < 55 TeV which corresponds to 
m Wi {m~ g ) < 400(1200) GeV! 

• Although we have not shown this here, we have checked that the same sign dilepton 
channel does not yield a better reach than the 4£ channel just discussed. Typically, 
cross sections in this channel are about 10-25% of the total dilepton cross section in 
Fig. |TT| whereas SM backgrounds (with just the basic cuts and triggers) are in the 
several fb range. 



D. Model Line D: The Higgsino NLSP scenario 

Within the GMSB framework described above, the value of that we obtain tends 
to be considerably larger than Mi and M 2 , so that the lightest neutralino is dominantly 
a gaugino (more specifically a bino, since Mi ~ |M 2 ). This is indeed the case for the 
three model lines examined above. Motivated by the fact that the phenomenology is very 
sensitive to the nature of the NLSP, we have examined a non-minimal scenario where we 
alter the ratio of \/j,\/Mi by hand from its model value, and fix a small value of \fi\ so that 
the NLSP is mainly Higgsino-like. We do not attempt to construct a theoretical framework 
which realizes such a scenario, but only mention that the additional interactions needed to 
generate \x and also the 5-parameter in this framework could conceivably alter the relation 
between \x and the gaugino masses. With this in mind, we use ISAJET to simulate a light 
Higgsino scenario where we take = 2, tan/5 = 3, M/A = 3, C grav = 1 but fix fi = — |Mi 
rather than the value obtained from radiative electroweak symmetry breaking. In practice, 
we do so by using the weak scale parameters obtained using the GMSB model in ISAJET as 
input parameters for the MSSM model, except that we use \i = — \M\ at this juncture. For 
this 'small //' scenario we expect that the two lightest neutralino and the lighter chargino 
will be Higgsino-like and close in mass, while the heavier charginos and neutralinos will be 



gaugino-like. The resulting spectrum is shown in Fig. |T3|a. Indeed we see that Z\ is the 
NLSP over the entire parameter range shown, and that Z 2 and W\ are generally only 20-30 
GeV heavier. As a result, the fermions from W\ and Z 2 decays to Z\ will be rather soft. 
Slepton masses are essentially family-independent because tan/3 is small. The lighter Higgs 
boson mass is just above 100 GeV, independent of A. 

Sparticle production at the Tevatron is dominated by the production of these Higgsino- 
like charginos and neutralinos as can be seen in Fig. [13]&. An important difference between 
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this case and chargino and neutralino production in the model lines examined above is that 
W\Z\ and Z\Z 2 production is also substantial. For a fixed chargino mass, however, the 
sparticle production cross section is somewhat smaller for model line C than it is for the 
other model lines. 

We have already noted that fermions from the decays W\ — > f f'Z\ and Z 2 — > f fZ\ are 
generally expected to be soft so that signatures for ZiZj or ZiWj production will closely 
resemble those for ZiZ\ production. In other words, sparticle signatures in such a scenario 



will be mainly determined by the Z 1 decay pattern shown in Fig. [L4|. For small values of 
A, Z\ — > G7. As A increases, the decays Z\ — > GZ and Z\ — > Gh become kinematically 
accessible, and the branching fraction for the photon decay becomes unimportant, while the 
decay to the Higgs boson becomes dominant. This is in sharp contrast to the gaugino NLSP 
case where the decay to the Higgs scalar is strongly suppressed. 

For small values of A (where Z\ mainly decays to via Z\ — > G7), the strategy for 
extracting the SUSY signal is as for model line A; i.e. to look for inclusive 27+ $t events. 
If we adopt the conservative background level of 1 fb as in this study, a '5cr' reach is obtained 
at the MI (TeV33) provided the signal cross section exceeds 3.5 fb (1 fb). These levels are 
shown as the horizontal dashed lines in Fig. [15| while the corresponding signal is shown by 
the curve labelled (7(77). We see that at MI (TeV33) experiments should be able to probe 
A values out to about 80 TeV (90 TeV) corresponding to ~ 165 GeV (180 GeV) via a 
search for 77+ $t events. 

For larger values of A, the NLSP dominantly decays via Z\ — > Gh and the di-photon 
signal drops sharply. In this case, since h mainly decays via h — > bb, the SUSY signal, which 
can contain up to four 6-quarks, will be characterized by multiple tagged 6-jet plus J^t events, 
which may also contain other jets, leptons and possibly photons (if one of the NLSPs decays 
via the photon mode). The dominant SM background to multi-6 events presumably comes 
from ti production and is shown in Table [IV], where we have also shown the signal cross 
section for A = 100 TeV. For events with one or two tagged 6-jets, the ti backgrounds come 
when the 6s from t decay are tagged; i.e. the rate for events where other jets are mis-tagged 
as 6's is just a few percent. This is also true for signal events. On the other hand, in the 
36 channel at least one of the tagged 6s in the ti background has to come from a c or light 
quark or gluon jet that is misidentified as a 6-jet, or from an additional 6 produced by QCD 
radiation. This is not, however, the case for signal events which contain up to four 6-jets. In 
each of the last two columns of Table |V| where we show the top background and the SUSY 
signal, we present two numbers: the first of these is the cross section when all the tagged 
jets come from real 6's, while the second number in parenthesis is the cross section including 
c and light quark or gluon jets that are mistagged as 6. Indeed we see that the bulk of the 
36 background is reducible and comes from mistagging jets, whereas the signal is essentially 
all from real 6s. 

It is clear from Table [IV] that the best signal to background ratio is obtained for events 
with > 36-jets. Our detailed analysis shows that although the signal cross section is rather 
small, 36-channel with a lepton veto (since top events with large $t typically contain leptons) 
offers the best hope for identifying the signal above SM backgrounds. We see that the signal 
is of similar magnitude as the background for A = 100 TeV, a point beyond the reach via 
the 77 channel. To further enhance the signal relative to the background we impose the 
additional cuts, 
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• #t > 60 GeV, and 

• 60 GeV < rribb < 140 GeV for at least one pair of tagged 6-jets in the event. 

The first of these reduces the signal from 2.7 fb to 2.1 fb while the background is cut by 
more than half to 3.4 fb. The mass cut was motivated by the fact that in these models, at 
least one pair of tagged 6's comes via h — > bb decay, with ~ 100 GeV, while the 6's from 
top decay form a continuum. We found, however, that this cut leads to only a marginal 
improvement in the statistical significance and the signal to background ratio. We traced 
this to the fact that, because of top event kinematics, one 6-pair is likely to fall in the 'Higgs 
mass window'. Reducing this window to 100 ±20 GeV leads to a slightly improved S/B but 
leads to too much loss of signal to improve the significance. 

The signal cross section via the 36 channel after the basic cuts as well as the additional 



$t and rribb cuts introduced above is shown by the solid curve labelled h — > bb in Fig. |T5. 
For small values of A, the signal is small because of the reduction in the branching fraction 
for Zi — > hG decay. The corresponding dashed lines shows the minimum cross section for a 
signal to be observable at the 5a level at the MI and at TeV33. We see that, at the MI, there 
will be no observable signal in this channel. In fact, even for the A value corresponding to the 
largest signal the statistical significance is barely 2a, so that a non-observation of a signal 
will not even allow exclusion of this model line at the 95% CL. With 25 fb~ l of integrated 
luminosity, the signal exceeds the 5a level for 82 TeV < A < 105 TeV (corresponding to 
rriyy ~ m z 1 ~ m z 2 ~ ^® GeV), and somewhat extends the reach obtained via the 77 
channel. Furthermore there appears to be no window between the upper limit of the 77 
channel and the lower limit of the 36-channel. A few points are worth noting. 

1. Since the background dominantly comes from events where a c or light quark or gluon 
jet is mis-tagged as a 6-jet, the reach via the 36 channel is very sensitive to our as- 
sumptions about the b mis-tag rate. Indeed, if the mis-tag rate is twice as big as we 
have assumed, there will be no reach in this channel even at TeV33. 

2. The 36 signal starts to become observable for A 80 TeV where the branching fraction 
for the decay Z\ — > hG becomes comparable to that for the decay Z\ — > 7G. The 
value of A for which the Higgs decay of the neutralino becomes dominant depends on 
rrih, which in turn is sensitive to tan/?. 

3. Although we have not shown it here, signals involving 6-jets together with additional 
photons or Z bosons identified via their leptonic decays have very small cross sections 
and appear unlikely to be detectable even for A ~ 100 TeV. 

Despite the fact that the top background alone is 50 to several hundred times larger than 



the SUSY signal in all relevant one and two tagged 6 plus multilepton channels in Table [TV], 
we have examined whether it was possible to separate the signal from the background. We 
focussed on the 26 + 0£ signal which has the best S/B ratio, and required in addition that 
$t > 60 GeV (which reduces the background by almost 50% with about a 20% loss of signal) 
and further 60 GeV < m&& < 140 GeV (which reduces the background by another factor 
of half with a loss of 25% of the signal) P5fl . We found, however, that the signal is below 
the 5a level over essentially the entire range of A even at TeV33: only for A = 80 ± 5 TeV 
does the signal cross section exceed this 5a level of 7.7 fb. Moreover the S/B ratio never 
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exceeds about 15% which falls below our detectability criterion S/B > 20%. We found that 
while it is possible to improve the S/ B ratio via additional cuts, these typically degrade the 
statistical significance of the signal. We thus conclude that for model line D, there will be 
no observable signal in the 2b channel even at TeV33. 

Before closing this discussion, it seems worth noting that we should interpret the reach 
in Fig. [15] with some care, because unlike in the study of model lines A, B and C, we do not 
really have a well- motivated underlying theory (that gives a Higgsino NLSP). We realized 
this by arbitrarily taking fi = — |Ml. The NLSP decay pattern, and hence the precise value 
of the reach, would depend on this choice which should be regarded as illustrative. In general, 
however, for the Higgsino NLSP model line, the coupling of the NLSP to Higgs bosons is 
substantial so that the branching fraction for the decay Z\ — > hG becomes large when this 
decay is not kinematically suppressed. For small values of A, such that the NLSP can only 
decay via Z\ — > jG, SUSY signals should be readily observable in the 77+ ~$t channel; once 
the NLSP decay to h begins to dominate, the cross section for diphoton events becomes 
unobservably small, and the SUSY signal mainly manifests itself as multiple b events which 
have large backgrounds from ti production. The most promising way to search for SUSY 
then seems to be via ft? events with > 3 tagged 6-jets but for a search in this channel an 
integrated luminosity of 25 fb^ 1 appears essential. A signal that extends the reach beyond 
that in the 77 channel is possible provided experiments are able to reduce the background 
from mis-tagged charm (light quark or gluon) jets to below 5% (0.2%). 

IV. SUMMARY AND CONCLUDING REMARKS 

The GMSB framework provides a phenomenologically viable alternative to the mSUGRA 
model. The novel feature of this framework is that SUSY breaking may be a low energy 
phenomenon. In this case, the gravitino is by far the lightest sparticle, and the NLSP decays 
within the detector into a gravitino and SM sparticles. Sparticle signals, and hence the reach 
of experimental facilities, are then sensitive to the identity of the NLSP. 

In this paper, we have examined signals for supersymmetric particle production at the 
Tevatron, and evaluated the SUSY reach of experiments at the MI or at the proposed TeV33 
within the GMSB framework. In our study, we consider four different model lines, each of 
which lead to qualitatively different experimental signatures. We use the event generator to 
simulate experimental conditions at the Tevatron, and for each model line, we have identified 
additional cuts that serve to enhance the SUSY signal over SM backgrounds. We assume 
the NLSP decay is prompt. This is a conservative assumption in that we do not make use 
of a displaced vertex to enhance the signal over SM background. 

For the first of these model lines, labelled A, the NLSP is mainly a hypercharge gaugino 
and dominantly decays via Z\ — > Gj, so that SUSY will lead to extremely striking events 
with multiple jets with hard leptons and large $t and up to two hard, isolated photons, 
with cross sections (after all cuts) shown in Fig. |j. The physics background to the 77 event 
topologies is very small, and detector-dependent instrumental backgrounds (such as from jets 
being mis-identified as photons or leptons, or large mismeasurement of transverse energies) 
dominate [|21j] . Even with a conservative estimate of 1 fb for the background cross section, 
experiments at the MI (TeV33) should be able to probe values of the model parameter A 
out to 110 TeV (130 TeV). If we optimistically assume that this background can be reduced 
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to a negligible level by hardening the cuts on the photons and we find a reach as high as 
A ~ 118 TeV (corresponding to a gluino of 950 GeV) at the MI and of 145 TeV at TeV33. 
For this model line, the ££77+ $t signal from slepton pair production may be observable 
at TeV33 even if sleptons are as heavy as 180 GeV. 

In model line B, the lighter stau is the NLSP, and heavier sparticles cascade decayed 
down to the stau, which then decays via f\ — > tG. The presence of isolated tau leptons 0, 
in addition to jets and other leptons is the hallmark of such a scenario. We found, however, 
that in some channels, the background from misidentification of QCD jets as r completely 
swamp the physics backgrounds, making it very difficult to detect the signal (see e.g. the 
Clrll channel in Table |) in this channel. We conclude that unless r misidentification can 
be greatly reduced from what we have assumed, SUSY will only be detectable via channels 
with at least three leptons (e, \x and r) for which the cross sections are individually small, 
and then, only by summing the signal in several channels. Even here, backgrounds from mis- 
identified taus are significant. Our assessment of the reach is shown in Fig. [| We see that at 
the MI the clean multilepton channel offers the best reach, out to A = 42 GeV (corresponding 
to = 192 GeV and squarks and gluinos as heavy as 700 GeV), while at TeV33, the 
reach may be extended out to A values somewhat beyond 50 TeV (trig ~ 800 GeV). 

In model line C, the lighter stau is again the NLSP but the lighter selectron and smuon 
are essentially degenerate with it, so that e.\ and [i\ cannot decay into a tau; i.e. these decay 
into a gravitino and a corresponding lepton. Since cascade decays of sparticles are equally 
likely to terminate in each flavour of slepton, we expect that this model line will lead to very 
large multiplicities of e, /1 and r in SUSY events. Indeed we found that the optimal strategy 
in this case was to search for events with n e + > 4 or n e + + n T > 4 with n e + n M > 2. 
The reach is shown in Fig. 12. We see that at the MI, there should be observable signals 
out to A = 45 TeV while at TeV33 A values as high as 55 TeV should be observable. These 
correspond to a gluino (chargino) mass of 1000 (320) GeV and 1200 (400) GeV, respectively! 

Finally, we have examined an unorthodox model line where, by hand, we adjust the 
value of fi to be smaller than the value of the hypercharge gaugino mass Mi. This leads 
to an NLSP which is dominantly a Higgsino. Furthermore, ~ m z 2 ~ m Zi so ^ na ^ 

the fermions from W\ and Z2 decays to Zi are soft, and SUSY event topologies are largely 
determined by the decay pattern of Z\. For small values of A for which the NLSP can only 
decay via Zi —>■ 7G, the signal is readily observable in the 77+ fix channel. However, 
once the NLSP decay to h begins to dominate, SUSY mainly manifests itself as multiple b 
events which have large backgrounds from tt production. The most promising way to search 
for SUSY then seems to be via fir events with > 3 tagged 6-jets and zero leptons. For a 
search in this channel, an integrated luminosity of ~25 fb^ 1 appears essential. A signal that 
extends the reach beyond that in the 77 channel is possible provided experiments are able 
to reduce the background from mis-tagged charm (light quark or gluon) jets to below 5% 
(0.2%). The reach for model line D is shown in Fig. [15], but it should be kept in mind that 
the details of this figure will be sensitive to our assumption about fi. 

To conclude, in GMSB models signals for SUSY events will be quantitatively and quali- 
tatively different from those in the mSUGRA framework. This is primarily because a neu- 
tralino NLSP decays into a photon, a Z boson or a Higgs boson and a gravitino, or sparticles 
cascade decay to a slepton NLSP which decays to leptons and a gravitino: these additional 
bosons (or their visible decay products), or leptons from slepton NLSP decays, often pro- 
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vide an additional handle which may be used to enhance the signal over SM background. 
Although we have not performed an exhaustive parameter scan, for the model lines that we 
studied we found that the SUSY reach (as measured in terms of the mass of the dominantly 
produced sparticles) is at least as big, and frequently larger than in the mSUGRA frame- 
work. For some cases, this conclusion depend on the capability of experiments to identify r 
leptons and 6-quarks with moderately high efficiency and purity. In view of the diversity of 
signals that appear possible for just this one class of models, we encourage our experimental 
colleagues to be in readiness for tagging third generation particles as they embark on the 
search for new phenomena in Run II of the Tevatron. 
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TABLES 



TABLE I. SM background cross sections in fb for various clean multilepton topologies from 
W, Z — > tt, VV (V = W, Z) and tt production at a 2 TeV pp collider, together with signal cross 
sections for A = 40 TeV and A = 50 TeV for Model Line B described in the text. For each event 
topology, the first number denotes the cross section after the basic acceptance cuts and trigger 
requirements along with the Z veto and the Acj) cut discussed in the text. The second number 
is after the additional cut, PTvis{T\) > 40 GeV, for events at least one identified r. The entries 
labelled Total* are the sum of all the cross sections except those in the 1t1£ channel. The last two 
rows provide a measure of the statistical significance of the signal. 
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TABLE II. SM background cross sections in fb for various jetty multilepton topologies from 
W, Z — > tt, VV (V = W, Z) and tt production at a 2 TeV pp collider, together with signal cross 
sections for A = 40 TeV and A = 50 TeV for Model Line B described in the text. The cross sections 
are with all the cuts including the prvis cut on the hardest r. 
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TABLE III. SM background cross sections in fb for clean and jetty events with > 4 leptons 
from V V (V = W, Z) production at a 2 TeV pp collider, after the basic cuts and triggers described 
in the text. At least 2 leptons are required to be e or fj,. Backgrounds from Z, W and tt production 
are negligible. Also shown are corresponding signal cross sections for A = 30 TeV and A = 40 TeV 
for Model Line C described in the text. As before, the C (J) refers to clean and jetty events. 
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TABLE IV. The background cross section in fb for multiple tagged 6-jets plus lepton plus 
fjT events from tt production after basic cuts and triggers discussed in the text. Also shown are 
the corresponding SUSY signal cross sections for A = 100 TeV for model line D. The numbers 
in parenthesis for the 36-channel include events from charm or light quark jets jets faking a 6-jet. 
Whereas this fake rate dominates the background in the 36 channel, it is negligible in the lb and 
26 channels. 
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FIG. 1. The four regions of the A — tan/3 parameter plane discussed in Sec. I of the text. 
We fix M = 3A, and take [i to be positive. The heavy solid lines denote the boundaries be- 
tween these regions. The grey region is excluded because electroweak symmetry is not correctly 
broken. The shaded regions should be probed by experiments at LEP and are nearly excluded 
because mf < 76 GeV (vertical shading), Z\ < 84 GeV (horizontal shading) or rrih < 95 GeV (or 
rriA < 83 GeV) (diagonal shading). The dot-dashed contours are where the chargino mass is 95, 
200 or 350 GeV, while the dotted line is the contour of = 90 GeV. 
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FIG. 2. a) Relevant sparticle masses, and b) production cross sections for the main sparticle 
production reactions at a 2 TeV pp collider versus the parameter A for the bino NLSP model line 
A. In frame b) the dot-dashed line represents the chargino pair production, while the dotted line 
denotes that for W1Z2 production. Also shown on the upper axis is the mass of the NLSP. 
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FIG. 3. a) The transverse energy distribution for the softer photon , and b) the $t distribution 
for the inclusive 77+ $t for model line A. 
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FIG. 4. Topological cross sections for inclusive a) O7, b) 17, and c) 77 plus fix signals after 
the cuts and triggers discussed in the text for model line A. The solid lines denote cross sections 
for events with at least one jet, while the dashed lines denote cross sections for events with no jets. 
The numbers on the lines denote the lepton multiplicity. The heavy solid line denotes the total 
signal cross section after all the cuts. 
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FIG. 5. Signal cross sections for clean ^£77+ $t events (£ = e, fx) from all SUSY sources (solid), 
and from just e.\ or \x\ pair production (dashed) versus A for model line A. The upper scale denotes 
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FIG. 6. The same as Fig. ^, but for the stau NLSP model line B. In frame a) the solid line 
denotes the mass of f\ while the dot-dashed line denotes the mass of the lighter selectron or smuon 
which are very nearly degenerate. 
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FIG. 7. The branching ratios for various decays of a) Z2, and b) l\ for model line B. The upper 
scale shows the mass of the parent sparticle. The relevant decay patterns of other sparticles are 
described in the text. 
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FIG. 8. Signal cross sections after all cuts versus A for the total signal in the clean (dashed 
line) and jetty (solid line) channels shown in Table | and Table O, respectively, for model line B. 
The corresponding horizontal lines denote the minimum cross section for a 5a signal, both at the 
MI and at TeV33. 
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FIG. 9. The same as Fig. |, but for the co-NLSP model line C. 
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FIG. 10. Branching ratios for decays of a) W\ and b) Z2 for model line C. Decay patterns of 
other sparticles are discussed in the text. 
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FIG. 11. Lepton multiplicity distributions for SUSY events after all cuts with A = 30 TeV, and 
A = 40 TeV for model line C. We show distributions of both n e + and n e + n^ + n T , where n T 
is the multiplicity of tagged r leptons in the event. 
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FIG. 12. Signal cross sections after all cuts for SUSY events with n e + n M + n T > 4 (solid) and 
n e + > 4 (dashed) for model line C. The corresponding horizontal lines denote the minimum 
cross section for a 5a signal (with a minimum of five signal events) at the MI and at TeV33. 
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FIG. 13. The same as Fig. ^, but for the Higgsino NLSP model line D. In frame a) the dotted 
line denotes the lightest neutralino, the solid line denotes the lighter chargino and the upper 
dot-dashed line denotes the second lightest neutralino. In frame b) the dashed, dot-dashed, solid 
and dotted lines denote cross sections for Z\Zi, W1Z2, W±Wi, W\Z\ production, respectively. 
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FIG. 14. Branching ratios for various decays of the neutralino NLSP versus A for model line 
D. The upper scale shows the neutralino mass. 
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FIG. 15. SUSY signal cross sections for the inclusive 77+ $t events (labelled (7(77), and for 
events with > 3 tagged 6-jets (labelled a(h — ► bb) after all cuts described in the text for model 
line D. The dashed horizontal lines denote the minimum cross section for the signal to be observable 
at the 5a level at the MI and at TeV33. 
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